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We present measurements of Ul and (j production at mid-rapidity from Au+Au collisions at 
nucleon-nucleon center-of-mass energies ^snn = 7.7, 11.5, 19.6, 27 and 39 GeV by the STAR 
experiment at the Relativistic Heavy Ion Collider (RHIC). Motivated by the coalescence formation 
mechanism for these strange hadrons, we study the ratios of N{Q~ -|- LI )/{2N{<j))). These ratios 
as a function of transverse momentum (pr) fall on a consistent trend at high collision energies, but 
start to show deviations in peripheral collisions at ^snn = 19.6, 27 and 39 GeV, and in central 
collisions at 11.5 GeV in the intermediate pr region of 2.4—3.6 GeV/c. We further evaluate em¬ 
pirically the strange quark pr distributions at hadronization by studying the Ll/<j) ratios scaled by 
the number of constituent quarks. The NCQ-scaled Ll/cj ratios show a suppression of strange quark 
production in central collisions at 11.5 GeV compared to ^snn > 19.6 GeV. The shapes of the 
presumably thermal strange quark distributions in 0-60% most central collisions at 7.7 GeV show 
significant deviations from those in 0-10% most central collisions at higher energies. These features 
suggest that there is likely a change of the underlying strange quark dynamics in the transition from 
quark-matter to hadronic matter at collision energies below 19.6 GeV. 

PACS numbers: 25.75.Dw, 25.75.Nq 


Lattice quantum chromodynamics (QCD) calculations 
suggest that, at high temperature and low baryon chem¬ 


ical potential ( rb ), the transition from the Quark Gluon 
Plasma (QGP) to the state of a hadron gas is smooth 


3 


and continuous (cross-over transition) [ij . At lower tem¬ 
peratures and high /tb, theoretical calculations predict a 
first order phase transition which may end at a crit¬ 
ical point (3|. The mapping of the QCD phase diagram 
has been a subject of intensive theoretical and experi¬ 
mental activities in the past decades. In central Pb-Pb 
collisions at Super Proton Synchrotron (SPS), the en¬ 
hanced production of fl at i^snn = 8-8 and 17.3 GeV 
and cj) mesons at y'snn = 6.3—17.3 GeV Q com¬ 
pared to TT mesons has been considered as a QGP sig¬ 
nal Multi-strange hadrons such as n(sss) hyperons 
and (j}{ss) mesons are importa nt p robes for the search 
of the QCD phase boundaries [lO, [Hi. The fl hyper¬ 
ons and (j) mesons are expected to have relatively small 
hadronic interaction cross sections [13, [il. Therefore, 
they can carry the information directly from the chem¬ 
ical freeze-out stage with little or no distortion due to 
hadronic rescattering. In addition, the measured D and 
(j) yields suffer minimal distortion from decay feed-down. 
As a result, the production of the D and (j) particles of¬ 
fers a unique advantage in probing the transition from 
partonic to hadronic dynamics. 

In heavy ion collisions at the top RHIC energy of .^snn 
= 200 GeV, model calculations |14l - [lq and experimen¬ 
tal data suggest that particles at intermediate pt are 
formed via the coalescence of low p^ quarks from the 
bulk partonic matter and/or fragmented hard partons. 
Experimentally, baryon to meson ratios have been found 
to be large compared to those from elementary colli¬ 
sions [l3, . The measured elliptic flow V 2 has been 

found to scale with the number of constituent quarks 
(NGQ) for both baryons and mesons [l^l in Au-|-Au col¬ 
lisions at the top RHIG energy. In order to explain 
these observations, coalescence model calculations re¬ 
quire the development of collectivity among constituent 
quarks during the partonic phase. This partonic collec¬ 
tivity has been considered as an important evidence for 
the formation of deconfined QCD matter with partonic 
degrees of freedom in Au-|-Au collisions at the highest 
RHIC energy [I3,[ii-l23. 

In order to map out the phase diagram of the QCD 
matter, a Beam Energy Scan (BES) program has been 
initiated at RHIC with Au-|-Au collisions at y^SNN = 
7.7 — 39 GeV [l^. These collisions allow us to reach a 
broad range of temperature and pB in tbe QCD phase 
diagram |25| and search for a possible beam energy region 
where the underlying dynamics are different from those 
of partonic matter observed in Au-|-Au collisions at the 
top RHIC energy. 

The STAR experiment, which has a large acceptance 
detector system [^, has collected Au-|-Au collision data 
at y^SNN = 7.7, II.5 and 39 GeV in 2010, and at 
y^SNN = 19.6 and 27 GeV in 2011. Compared to previ¬ 
ous measurements at SPS [H, RHIC offers wide range 
of beam energies. In addition, the STAR detector pro¬ 
vides uniform acceptance over different energies and ex¬ 
tensive reach to the intermediate pt range for both D 
and (j) mesons for different centralities. In this letter. 


we present the first STAR measurements of mid-rapidity 
(|y| < 0.5) D and (j) production with a broad pr coverage 
for various collision centrality bins at selected energies 
for the BES. With guidance from coalescence formation 
mechanism for primordial particles such as D and cj), we 
examine features in particle production and explore pos¬ 
sible changes from bulk partonic coalescence to hadron 
dominated dynamics as the colliding energy decreases. 

A minimum bias trigger was used to record Au-fAu 
collision events in this analysis. The trigger was defined 
using a coincidence of signals from either the Zero Degree 
Calorimeters. Vertex Position Detectors or Beam-Beam 
Counters [23, [11|. STAR’S Time Projection Chamber 
(TPC) was used for tracking of charged particles and 
particle identification. In the offline data analysis, we re¬ 
quire the radial position of the reconstructed primary 
vertex to be within 2 cm of the beam axis to suppress 
events from collisions with the beam pipe (radius of 3.95 
cm). To ensure nearly uniform detector acceptance, the 
analyzed events were required to have a primary Z vertex 
(along beam direction) within ±70 cm from the center of 
the TPC for ^snn = 7.7, 19.6, and 27 GeV collisions and 
±50 cm, ±40 cm for ^snn = 11-5 and 39 GeV, respec¬ 
tively. After the event selection, we obtain approximately 
4, 12, 36, 70 and 130 million Au±Au minimum bias trig¬ 
gered events at ^snn = 7.7, 11.5, 19.6, 27 and 39 GeV, 
respectively. The collision centrality was determined by 
comparing the uncorrected charged hadron multiplicity 
measured from the TPC at mid-rapidity (|? 7 | < 0.5) with 
Monte Carlo Glauber simulations (TtI. [^. 

The multi-strange hadron signals and raw yields were 
obtained from the invariant mass distributions recon¬ 
structed through their hadronic decay channels: (j) —>■ 
K+ + K- and D-(O^) A(A) ± R-{!<+). The de¬ 
cay daughters A(A) were reconstructed through A(A) —>■ 
p{p) ± 7r“(7r“*'). Charged hadrons {t:^,K^,p,p) were 
identified by their specific energy loss {dE/dx) in the 
TPC gas [^. The combinatorial background of the 
weakly decaying particles A(A) and D“(D''') was reduced 
by geometrical cuts on their decay topology [H, [s^ . The 
fl~ ) combinatorial background was estimated by ro¬ 
tating K~{K^) tracks at 5 different angles from 7r/3 to 
57r/3 and normalizing the invariant mass distribution to 
the mass window of (1.625 GeV/c^, 1.655 GeV/cP) and 
(1.69 GeV/c^, 1.72 GeV/c^). The raw yields 

were extracted by counting the signals within a mass 
window from 1.660 to 1.685 GeV/c^ after subtracting 
the rotational background. The K'^K~ combinatorial 
background in (j) meson reconstruction was subtracted 
with the mixed event technique [T3 . l3lj| . The (j) meson 
raw yields were determined by a Breit-Wigner ± poly¬ 
nomial function (up to second order) fit to the mixed- 
event-background-subtracted K~ invariant mass dis¬ 
tribution [io, Hlj. 

Figures 1 and 2 show the pr spectra of (j), and 
at mid-rapidity (|?/| < 0.5) for different centralities from 
Au±Au collisions at ^snn = 7.7 — 39 GeV. The spec- 
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ripheral collisions by 2%-3% due to greater combinatorial 
backgrounds. For px < 0.8 GeV/c in central collisions, 
the uncertainty of (j) meson raw yield extraction is dom¬ 
inant. However, for px > 1.6 GeV/c the main source 
of systematic error is the differences in track selection 
cuts. For the H invariant yields, the relative systematic 
errors vary from ^ 5% to 20%, and are dominated by the 
signal extraction methods. Due to the higher combina¬ 
torial background in px 1.2 GeV/c and low statistics 
at Px ^ 2.8 GeV/c, the systematic errors are found to 
be larger in the corresponding px windows. The system¬ 
atic uncertainties have a weak centrality dependence and 
their energy dependences for Q and (p particles are simi¬ 
lar. The systematic errors of invariant yields of (j) and D 
are shown as green bands in Figs. 1 and 2 for each px 
bin. 


FIG. 1: (Color online) Mid-rapidity (|i/| < 0.5) (j) meson px 
spectra from Au-|-Au collisions at different centralities and 
energies (-^snn = 7.7 - 39 GeV). The green bands represent 
systematic errors. The dashed curves represent fits to the 
experimental data with a Levy function (ICll . 


■a 

H 

a 


1 

10'^ 

10-^ 

10 -‘ 

10 -* 

10-‘» 

1 


z 


10 " 

10 * 

lo” 

10 '“ 


Au+Au 7.7 GeV ' 

Au+Au h.5 GeV ■ 

' Au+Au 19.6 GeV ' 

■ 

: n 

*++ 

- 

oi + 

■ * 

- 


+ 0-10% 


* 0-60% 


r.- 

▼ 10%-60%/10 


Au+Au 27 GeV 

Au+Au 39 GeV 

Solid: a' 

4 - 

- izHit 

^ * 

- * 

* 

- ®8*l**«:l* 

■ 

Open:IT^/2 

19.6 - 39 GeV 
+ 0-10% 

• 10%-20%/10 
■ 20%-40%/10^ 

: ■+■ 

* 40%-60%/10’ 

* 60%-80%/10-‘ 


012345 12345 12345 

Transverse momentum (GeV/c) 


FIG. 2: (Color online) Mid-rapidity (|t/| < 0.5) D“(fl^) px 
spectra from Au-|-Au collisions at different centralities and 
energies (.^snn = 7.7 - 39 GeV). The open symbols represent 
and solid symbols represent The green bands denote 
systematic errors. 


tra were corrected for reconstruction efficiency and ge¬ 
ometrical acceptance. The systematic errors mainly 
come from two sources: the different signal extraction 
techniques, and the reconstruction efficiency corrections. 
They were studied as a function of px and were obtained 
by exploring the dependence of invariant yields on var¬ 
ious raw yield extraction techniques including different 
fit/counting ranges and different fit functions, and on 
different combinations of analysis cuts. For the </> me¬ 
son, relative systematic errors of invariant yields vary 
from 10%-I6% at = 11.5 - 39 GeV to 17%-21% at 
^snn = 7-7 GeV. The systematic errors in 0-10% central 
collisions are generally larger than those in 60%-80% pe¬ 



FIG. 3: (Color online) The baryon-to-meson ratio, N -|- 
Q )/{2N{(j})), as a function of px in mid-rapidity (|i/| < 0.5) 
from Au+Au collisions at y'SNN = 7.7 — 200 GeV. The green 
bands denote systematical errors. The solid and dashed lines 
represent recombination model calculations for central colli¬ 
sions at .^snn = 200 GeV [l^ with total and thermal strange 
quark contributions, respectively. 

We present baryon-to-meson ratios, N{il~ + 
fl^)/{2N{(j))), as a function of px from Au+Au 
collisions for various beam energies from y^snn = 7.7 to 
200 GeV in Fig. [3] and for various collision centralities 
in Fig. 01 respectively. Data from 200 GeV Au+Au 
collisions are from previously published STAR results 
[a. Coalescence or recombination models [iM3 have 
been used to describe particle productions in nucleus- 
nucleus collisions at RHIC. In particular, a model 
calculation by Hwa and Yang for Au+Au collisions at 
y^NN = 200 GeV [l^ predicted that most of the D 
and (j) yields up to the intermediate px region are from 
coalescence/recombination of thermal strange quarks. 
The straight dotted line assumed that these thermal 
strange quarks have exponential px distributions. Devi¬ 
ations from the straight line at high px were attributed 
to recombination with strange quarks from high px 
showers. Deviations from the theory calculation at low 
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Pt could indicate that thermal strange quarks may not 
have an exponential distribution. Possibly, other particle 
production dynamics may also contribute deviations 
from coalescence model calculations. 

In Fig. 13] the measured N{Q~ + ratios 

from central Au+Au collisions at y^SNN = 19.6, 27 and 
39 GeV follow closely the ratio from 200 GeV and are 
consistent with a picture of coalescence/recombination 
dynamics over a broad px range of 1 — 4 GeV/c. The ra¬ 
tios at 11.5 GeV seem to deviate from the trend observed 
at higher beam energies. In particular, the ratios at 11.5 
GeV appear to turn down around pr of 2 GeV/c while 
those at higher beam energies such as 39 and 200 GeV 
peak at px of 3 GeV/c or above. The collision central¬ 
ity dependence of the N{n~ +0.'^)/{2N ratios in Fig. 
|4Ua)-(d) shows distinct differences between the 40%—60% 
centrality bin and the other centrality bins for Au-t-Au 
collisions at 19.6 and 27 GeV. Furthermore, the ratios 
from the peripheral collisions of 40% —60% at 27 GeV are 
similar in magnitude to the ratios from collisions at 11.5 
GeV. Because the 17 and 4> particles have small hadronic 
rescattering cross sections [3^ , the change in these ratios 
is likely to originate from the partonic phase. The de¬ 
crease in the N{n~ + )/{2N{(j))) ratios from central 

collisions at 11.5 GeV compared to those at 19.6 GeV 
or above may indicate a significant change in the hadron 
formation dynamics and/or in strange quark px distribu¬ 
tions at the lower energy. Such a change may arise from 
a transition from hadronic to partonic dynamics with in¬ 
creasing beam energy. The turn-over in the ratios from 
Au-|-Au collisions below 11.5 GeV beam energy is un¬ 
likely due to contributions of high px shower partons as 
suggested by model calculation from Hwa and Yang [l^ 
because of relatively low px particles involved. 


We illustrate qualitatively the change in the underly¬ 
ing bulk strange quark px distribution by following the 
procedure developed in Ref. [l^. We assume that the 
il baryons are formed from coalescence of three strange 
quarks of approximately equal momentum and the 4> 
mesons from two strange quarks. In the coalescence 
framework, the 17 baryon production probability is pro¬ 
portional to the local strange quark density, fs{px), and 
the (p meson is proportional to fs{Px)f's{Px): where fsifs) 
is the strange (anti-strange) quark px distribution at 
hadronization. Assuming that strange quarks and anti- 
strange quarks have a similar slope of px distribution, 
Ar(a--Hn+)| 0^3 

the NCQ-scaled ratio —- -—— could reflect the 

strange quark distribution at hadronization. We note 
that theoretical calculations with a more sophisticated 
recombination scheme have been developed by He et al. 
[ 3 ^ and the extracted strange quark distribution is sim¬ 
ilar to that from the present approach. 


Figure jS^a) shows the NGQ-scaled 


NiQ~+n+)\^n^^ 

2Af(0)l>_L 


ratios as a function of p^ = px/^q at mid-rapidity 
(|y| < 0.5) from central Au-I-Au collisions at y^snn = 



FIG. 4: (Color online) Centrality dependence of A(17 -|- 

17 )/{2N{(j))) ratios, as a function of px in mid-rapidity 
(|y| < 0.5) from Au-|-Au collisions at ^snn = 11.5, 19.6, 
27 and 39 GeV. The green bands denote systematical errors. 


11.5 — 200 GeV as well as 0-60% collisions at 7.7 GeV. 
Here iV(f7“ -|- 17 ) and denote invariant yields of 

(17“ -1-0) and (j>, respectively. Since the px bin widths 
used for the 0”(0~'") and (p meson spectra do not match, 
we use our Levy fit (see Fig. [T|) to interpolate the in¬ 


variant yield of 4> meson at desired px- The NCQ-scaled 


jv(n"+a+)|pn=3 


ratios at all energies can be fit with a 


Boltzmann distribution g (mr 

T(ms+T) ’ 

is the effective strange quark mass of 0.46 GeV/c^ from 


Ref. [ 1 ^, mx is the transverse mass {y/mj + Px), T is 
the slope parameter of the exponential function which 
may be related to the freeze-out temperature and radial 
expansion velocity of strange quarks [solj . Considering 
different yields ratios of s quark over s quark with col¬ 
lision energies, that is, fs{Px) = '''{^/sNN)fs{Px)^ where 


?’^(i/snn) = ^(17'*')/^ (17 ), we include a correction 


factor ps = (l-|-r^)/r’ in the Boltzmann distribution func¬ 
tion (based on the coalescence calculation [H) , and then 
A is proportional to strange quark rapidity density. 


The fitting parameters A and T, and Icr contours are 
shown in Fig. ISKb). Figs. |S](a)-(b) show that the de¬ 
rived strange quark distributions vary little in shape as 
a function of beam energy from 11.5 GeV to 200 GeV. 
The amplitude parameter A at 11.5 GeV, however, seems 
to be noticeably smaller than those data of 19.6 GeV 
or above. Based on coalescence model [l^, the smaller 
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0.002 0.004 0.006 0.008 

A oc dN/dy of s quark 



FIG. 5: (Color online) (a) NCQ-scaled N{n- +0.'^)/{2N{(!>)) 
ratios, as a function of in mid-rapidity (|t/| < 0.5) 

from Au-|-Au collisions at ^snn = 7.7 — 200 GeV. Here Uq is 
the number of constituent quarks of each hadron. The green 
bands denote systematic errors. Dashed lines are Boltzmann 
fits to data, (b) The fitting parameters A and T, and la 
contours (including statistic and systematic errors), (c) The 
ratios of /s(Pt) and scaled ^ function of p^. 

The scale factors are 394.4, 763.7, 742.8, 870.9, and 746.5 for 
^snn = 7.7, 11.5, 19.6, 27, and 39 GeV, respectively. 


strange quark local density at 11.5 GeV is probably re¬ 
sponsible for the smaller N{n~ + )/{2N{(/))) ratios as 
shown in Fig. 3, where the first two low px points at 11.5 
GeV are systematically lower than those at y^SNN > 19-6 
GeV. At 7.7 GeV, the slope parameter T is smaller than 
those data of 19.6 GeV or above, with a 1.8cr standard 
deviation from the 19.6 GeV result. We note that one 


possible reason for the deviation of T is the central¬ 
ity difference since the data at 7.7 GeV are for 0-60% 
while those at other energies are for central collisions. 
In the framework of the coalescence mechanism, our de¬ 
rived ratio distribution can be sensitive to both the lo¬ 
cal density and the px distribution of strange quarks. 
Our data of 19.6 GeV or above show little beam energy 
dependence suggesting strange quark equilibration may 
have been approximately achieved in those central colli¬ 
sions, possibly due to strange quark dynamics rather than 
hadronic processes [s^. The variation of the 11.5 GeV 
data may arise from the strangeness non-equilibration 
and the presence of a strangeness phase space suppres¬ 
sion factor (ys < 1) [s^. A possible transition in the 
collision dynamics and in the dominate degrees of free¬ 
dom (partonic versus hadronic) below 19.6 GeV needs 
further experimental investigation with more experimen¬ 
tal probes and with larger data samples [ssj j. 

Recently, the ALICE experiment reported an observa¬ 
tion of nearly flat ratios of proton to (/) as a function of px 
in central Pb-|-Pb collisions at .^snn = 2.76 TeV [s^. It 
was argued that the similarity in shapes of the px spectra 
indicates that the radial flow of these particles is mostly 
determined by the masses of these particles as hydrody- 
namical calculations predicted, instead of by the number 
of constituent quarks as expected from coalescence mod¬ 
els [1^. The ALICE proton to (p ratios as a function 
of Px also showed a very strong dependence on collision 
centrality. We note that such a strong dependence is an 
indication that the protons undergo considerable rescat¬ 
tering during the hadronic evolution as expected from 
the hybrid model calculation [s^. It is important to dis¬ 
entangle the hadronic rescattering contributions to the 
radial flow for ordinary hadrons in order to address the 
partonic flow prior to hadronization. We have used the 
n and (j) spectra to carry out another independent check 
on the coalescence picture: we can obtain the strange 
quark Px distribution from the quark number scaled 
to (j) ratios as in Eig. [U^a); and we can calculate another 
strange quark px distribution by dividing the (p with the 
previously obtained strange quark distribution. Fig.jSKc) 
shows these two strange quark distributions. The con¬ 
sistency of these distributions indicates that there is one 
unique strange quark distribution which can explain both 
n and p Px spectra, a necessary condition for coalescence 
model. 

In summary, STAR has measured the production of 
multi-strange hadrons fl and p at mid-rapidity from 
Au-I-Au collisions at y'snn = 7.7, 11.5, 19.6, 27 and 
39 GeV from the BES program at RHIG. The A(fl“ -|- 

12^)/{2N{p)) ratios at intermediate pt in peripheral col¬ 
lisions are found to be lower than those in central colli¬ 
sions at 19.6, 27 and 39 GeV. The ratios from 11.5 GeV 
central collisions are systematically lower than those from 
collisions at 19.6 GeV or above for px > 2.4 GeV/c. The 
NCQ-scaled ft/p ratios show a suppression of strange 
quark production in 11.5 GeV compared to y^SNN > 19.6 
GeV. The shapes of the presumably thermal strange 
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quark distributions in 0-60% most central collisions at 
7.7 GeV show significant deviations from those in 0-10% 
most central collisions at higher energies. These features 
suggest that there is likely a change in the underlying 
strange quark dynamics in the bulk QCD matter respon¬ 
sible for n and (j) production. Our measurements point 
to collision energies below 19.6 GeV for further inves¬ 
tigation of a possible transition from partonic dominant 
matter (i/snn > 19.6 GeV) to hadronic dominant matter 
(v«< 11-5 GeV). 
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